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ABSTRACT 

We explore the spatio-temporal evolution of solar flares by fitting a radial expansion model 
r{t) that consists of an exponentially growing acceleration phase, followed by a deceleration phase 
that is parameterized by the generalized diffusion function r{t) oc /^(t — ti)^/^, which includes the 
logistic growth limit (/3 = 0), sub-diffusion (/3 = — 1), classical diffusion (/3 = 1), super-diffusion 
(/3 = 1 — 2), and the linear expansion limit (/3 = 2). We analyze all M and X-class flares observed 
with GOES and AIA/SDO during the first two years of the SDO mission, amounting to 155 
events. We find that most flares operate in the sub-diffusive regime {/3 = 0.53 ± 0.27), which we 
interpret in terms of anisotropic chain reactions of intermittent magnetic reconnection episodes in 
a low plasma- /3 corona. We find a mean propagation speed of v = 15± 12 km s~^, with maximum 
speeds of Vmax = 80 ± 85 km s~^ per flare, which is substantially slower than the sonic speeds 
expected for thermal diffusion of flare plasmas. The diffusive characteristics established here (for 
the first time for solar flares) is consistent with the fractal-diffusive self-organized criticality (FD- 
SOC) model, which predicted diffusive transport merely based on cellular automaton simulations. 

Subject headings: Sun: Solar Flares — Statistics — magnetic fields 



1. INTRODUCTION 

Many physical processes can be characterized by their spatio-temporal evolution, which we simply de- 
fine here as the temporal evolution or time dependence of a spatial or geometric parameter, say x(t), where 
the geometric parameter x could be a spatial length scale, an area, a volume, a fractal dimension, or any 
combination of these. At the largest scales, for instance, the Big Bang theory describes the evolution of 
the universe by the gradual expansion of its radius r(t), which can be decelerating in a closed universe 
(Friedmann-Lemaitre- Robertson- Walker model). Recent supernova observations reveal an accelerating uni- 
verse at z ~ 0.5 (Riess et al. 1998; Perlmutter et al. 1999), and a deceleration that preceded the current 
epoch of cosmic acceleration (Riess et al. 2004). 

Other dynamic processes that predict a specific spatio-temporal evolution based on some physical model 
include Brownian motion, fractal Brownian motion, classical diffusion, sub-diffusion, super-diffusion. Levy 
flights, logistic growth, percolation, self-organized criticality avalanches, cellular automatons, non-extensive 
Tsallis entropy, complex networks, etc. The knowledge of the spatio-temporal evolution of a physical process 
is often directly related to scaling laws between various physical parameters, and thus spatio-temporal 
measurements play a decisive role in the derivation of scaling laws. For instance, particle conservation 
{nV = const) in an expanding gas, plasma, or universe predicts a reciprocal scaling between the particle 
density n{t) and the volume V{t), i.e., n{t) oc V{t)~^^ and consequently a scaling of n{t) oc r{t)~^ as a 
function of the radius r{t) in the case of homogeneous isotropic expansion, while it scales as n{t) oc r{t)~^ 
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for a fractal volume with Hausdorff dimension D < 3. The existence of a spatio-temporal evolution, such as 
n{t) oc r{t)~^ for adiabatic expansion, implies then also the prediction of a statistical correlation or scaling 
law n oc r~^, if the expansion speeds v = dr/dt of the sample have a limited range. 

Needless to say, that scaling laws obtained from solar data, where we have ample spatial resolution, 
are extremely useful for the interpretation of stellar data, where we have no spatial resolution at all and 
have to rely on scaling laws measured in solar or magnetospheric plasmas. Here we focus on the spatio- 
temporal evolution of solar flares, which is a completely unexplored topic, but bears important information 
on the underlying dynamic processes. There are very few statistical measurements of spatial scales L, areas 
A, and volumes V, or fractal dimensions D of solar flares, and virtually no statistical studies about the 
temporal evolution of these parameters, such as L{t)^ A{t)^V{t)^ D{t). A few statistical measurements of 
spatio-temporal parameters of solar flares (compiled in Aschwanden 1999) have been made from the S- 
054 soft X-ray imager onboard Skylab (Pallavicini et al. 1977), from Yohkoh soft X-ray images (Kano and 
Tsuneta 1995; Porter and Klimchuk 1995; Aschwanden et al. 1996; Metcalf and Fisher 1996; Reale et al. 1997; 
Shimizu 1997; Garcia 1998; Nagashima and Yokoyama 2006), from the Multispectral Solar Telescope Array 
(MTSA) rocket flight (Kankelborg et al. 1997), from the extreme- ultraviolet (EUV) imager SOHO/EIT 
(Berghmans et al. 1998; Krucker and Benz 2000), and from the EUV imager on TRACE (Aschwanden 
et al. 2000; Aschwanden and Parnell 2002; Aschwanden and Aschwanden 2008a, b). Most of these studies 
provide statistics on spatial length scales L, areas A, and durations T of flares (down to nanoflares), but 
there exists no study to our knowledge that provides statistics on the spatio-temporal evolution L(t) or A{t) 
of solar flares. 

In this paper we are going to analyze the spatio-temporal evolution of all large (GOES X- and M- 
class) flares observed during the first two years of the SDO mission, which were observed with high spatial 
resolution (0.6^^), high cadence (12 s), and in 7 coronal wavelengths filters that cover a wide temperature 
range (T ^ 0.5 — 16 MK). This is an ideal data set for such a study, because both spatial and temporal 
parameters can be measured with unprecedented accuracy and 100% time coverage. Setting a threshold of 
>M1.0 GOES class, we obtain a complete set of 155 flare events detected with both GOES and AIA, which 
makes it to a perfect representative statistical sample. The data analysis presented here is restricted to the 
335 A filter, which appears to be very suitable to capture the high-temperature component of the energy 
release and heating phase of these largest flares. We fit then various theoretical models to the measured time 
profiles of the fiare areas A{t) or mean radius r{t) oc ^(t)-^/^, such as: classical diffusion, r{t) cx t^/^; sub- 
diffusion r(t) oc t^/^ with /3 < 1; super-diffusion or Levy flights, r(t) oc t^/^ with (3 > 1; or logistic growth, 
where the time profile initially expands exponentially and then saturates at a finite level, a limit that is also 
called carrying capacity of limited resources in ecological models. The generic form of these spatio-temporal 
evolution functions r{t) are shown in Fig. 1. From the fits of the theoretical models to the observed flare 
data we aim then to gain physical insights into the underlying dynamic processes, which concern the spatial 
propagation (or chain reaction) of magnetic reconnection or other nonlinear energy dissipation processes. 
A particular statistical model that we test is the so-called fractal-diffusion self-organized criticality model, 
which predicts specific spatio-temporal evolutions and powerlaw distributions of the spatial and temporal 
parameters. 

The plan of this paper consists of a brief description of relevant theoretical models (Section 2), the 
statistical data analysis and forward- fitting to 155 flare events observed with AIA/SDO (Section 3), a 
discussion of the interpretation and consequences of the results (Section 4), conclusions (Section 5), and a 
generalization of the FD-SOC model (Appendix A). 
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2. THEORETICAL MODELS 

In the following we introduce some general theoretical models that describe the spatio-temporal evolution 
of physical systems, which have some universal validity for the dynamics of a large number of nonlinear and 
complex systems. The reason for their universal validity and applicability lies in the fact that they are 
formulated only in terms of the fundamental parameters of space and time, while individual applications 
generally involve more specific physical parameters (such as densities, temperatures, or magnetic fields in 
astrophysical applications) . 



2.1. Classical Diffusion and Brownian Motion 

Classical diffusion obeys the following differential equation, 

df{x,t) d^f{x,t) 

where f{x^t) is the spatio-temporal distribution function of particles and n is the diffusion coefficient. The 
diffusion equation, expressed here for a 1-dimensional (1-D) space coordinate x, can be generalized to 2-D 
or 3-D space, where x then signifies the distance r from the center position of the initial distribution. The 
solution of the diffusion equation is a Gaussian distribution function (Einstein 1905), 

with a time- dependent evolution of its second moment, 

(x^) = 2K.t , (3) 
which simply means that its radial size r{t) = \/~(^^ evolves proportionally to the square root of the time, 

r(t) oc , (4) 

since the start of the process from an initial ^-function. This spatio-temporal relationship has originally 
been applied to the molecular motion in a fluid (Brown 1828), and thus is also called Brownian motion, 
named after the Scottish botanist Robert Brown. Experimentally it can be demonstrated by watching the 
vibrations of a dust particle suspended in a fluid under a microscope. While the random walk of a single 
particle is stochastic and unpredictable at every time step, the average distance of the motion of an ensemble 
of particles or the radial expansion of the ensemble is well-characterized by the square-root relationship of 
Eq. (4), and is also called classical diffusion, or more popularly "the random walk of a drunkard". Numerous 
applications of classical diffusion can be found in nature, laboratory physics, and astrophysics, such as the 
expansion of released gases, aerosols in the atmosphere, dust clouds from vulcanos, plasma diffusion in fusion 
reactors, thermal diffusion in solar flare plasmas, or insect ecology in biophysics. For a recent textbook on 
applications of Brownian motion see, e.g., Earnshaw and Riley (2011). 
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2.2. Anomalous Diffusion, Sub-DifFusion, Super-Diffusion, Levy Flights 

Deviations from classical diffusion processes, also called anomalous diffusion, are often defined in terms 
of a non- linear dependence on time, characterized with a power law index /3 deviating from unity, 

{P < 1 (sub-diffusion) 
f3 = 1 (classical diffusion) (5) 
(3 > 1 (super-diffusion or Levy flight) 

We show the generic time evolution of a sub-diffusion process with f3 = 1/2 and a super-diffusion process 
with /3 = 3/2 in Fig. 1. Anomalous diffusion implies more complex properties of the diffusive medium than 
a homogeneous structure, which may include an inhomogeneous fluid or fractal properties of the diffusive 
medium. Anomalous diffusion was found in biology, such as active cellular transport, protein diffusion within 
cells, or diffusion through porous media (percolation phenomena). 

Related to super-diffusion processes are Levy flights, named after the French mathematician Paul Pierre 
Levy by Mandelbrot (1982), which are defined similarly to super-diffusion processes (with /3 > 1 in Eq. 5) 
and lead to heavy-tailed (powerlaw) distribution functions /(x,t), in excess of the Gaussian distribution 
functions (Eq. 2) obtained for classical diffusion. A popular example for Levy flights is the random walk of 
sharks or other ocean predators, which abandon Brownian motion by occasional large-range jumps, when 
they do not find sufficient food. Other applications of Levy walks involve diffusion on fractal structures, 
tracer diffusion in living polymers, turbulent rotating flows, electromagnetically driven flows, subrecoil laser 
cooling, chaos in a Josephson junction (electronic chip), or even the geographic travel of bank notes. For an 
overview see, e.g., Zumofen et al. (1999). 



2.3. Fractal-Diffusion in Self-Organized Criticality Systems 

Inhomogeneous media can be structured in complex patterns. One of the simplest complex patterns 
are fractal structures (Mandelbrot 1982), which can be characterized by a single parameter, such as the 
fractal dimension D. For instance, a cubic volume V with length L that has a space-filling Euclidean volume 
y = L^, can have a fractal substructure with a fractal volume V = L^, with D < 3 being the (fractal) 
Hausdorff dimension. 

Such fractal volumes were found to describe well the avalanches in nonlinear dissipative systems in the 
state of self-organized criticality (SOC), which led to the statistical fractal- diffusive model of a slowly-driven 
self-organized criticality system (Aschwanden 2012). The prototype of a SOC model is the cellular automaton 
model, which numerically simulates the spatio-temporal evolution of a SOC avalanche by thresholded next- 
neighbor interactions. The "spatial diffusion" of a SOC avalanche is defined by an iterative mathematical 
re-distribution rule between next neighbors in a grid (i, j. A:), such as, 

z{i,j, k) = z{i,j, k) -\- 1 initial input 

z{i, j, k) = z{i, j, k)-S if z{i, j, k) > 8, (6) 

z{i±lj ±l,k±l) = z{i±lj ±l,k±l) + 1 

where z denotes the number of energy quanta in each node that are re-distributed in the nonlinear energy 
dissipation process of a SOC avalanche (Bak et al. 1987). It was found that the instantaneous volume change 
dV{t)/dt of a SOC avalanche can be approximated with a fractal scaling (Aschwanden 2012), 

(T) 
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and evolves in space according to the classical diffusion law (Eq. 4) 

r{t) oc t^/' . (8) 

This scaling law has been inferred from the spatio-temporal evolution of a few cellular automaton avalanches 
(see Fig. 3 in Aschwanden 2012), which constitutes a strong prediction for any observation of SOC systems 
in nature. This prediction, however, has never been tested with observational data, such as for solar flares, 
which we set out to test in this study for the first time. 



2.4. Logistic Growth Model 

Instabilities represent a loss of equilibrium and show often a time evolution that consists of an initial ex- 
ponential growth phase (which we might call the "acceleration phase"), followed by quenching and saturation 
of the instability in the decay phase (which we might call "deceleration phase" ) . The total dissipated energy 
in such an instability is often approximately proportional to the unstable volume. We can then describe 
the exponential growth phase with a proportionality of the change of volume dV{t)/dt to the instantaneous 
volume F(t), 

where V is the growth rate, and tg = l/T denotes the e- folding growth time, which has the simple solution 
of an exponential function, 

y(t) = Foexp( — ) . (10) 

TG 

A more general approach to describe both the exponential growth phase together with the saturation phase 
is the logistic equation^ which has been widely used in ecologic applications, but has universal validity for 
nonlinear systems with limited free energy. The logistic equation is defined by a simple first-order differential 
equation, discovered by Pierre Francois Verhulst in 1845 (see textbooks on nonlinear dynamics, e.g.. May 
1974; Beltrami 1987, p.61; Jackson 1989, p.75; Aschwanden 2011, p.94). 



dV{t) _ V{t) 

dt TG 



1 v{t) 



(11) 



where Voo represents the asymptotic limit that is reached in the saturation phase. In ecologic applications, 
the asymptotic limit Voo is also called carrying capacity^ such as the total amount of available energy or 
ressources that can be sustained world-wide. From Eq. (11) we see that the logistic equation approaches the 
exponential growth equation (Eq. 9) for small times (when V{t) <C Voo), and the limit dV{t)/dt ^ and 
V{t) ^ Voo for large times. The explicit solution of Eq. (11) is, 

V{t) = , , , , (12) 

l + exp(-^) ^ ^ 

while the time derivative dV/dt (representing the energy dissipation rate dE/dt if the energy is proportional 
to the volume), 

dV(t) _ Vq^ exp(-^^) 



dt TG [l+exp(-^^)]2 • 

If we approximate the unstable volume with a spherical geometry, the spatial length scale or radius r{t) is 
related to the volume V{t) as, 

V{t) = jnr{t)' , = jnrl , (14) 
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which predicts the following spatio-temporal evolution r{t) 

3 



r{t) 



47r 



V{t) 



1/3 



1 + exp(- 



t-ti 

TG 



1-1/3 



(15) 



A graphical illustration of this logistic growth function is shown in Fig. 1, which appears to conform to a 
lower limit of sub- diffusion processes (for (3 0). The logistic growth model was applied to population 
growth, neural networks, tumor growth in medicine, autocatalytic reactions in chemistry, the statistical 
distribution of fermions in atomic physics, language change in social sciences, or diffusion of innovations in 
economics. In solar physics, the logistic-growth model was applied to model hard X-ray pulses (Aschwanden 
et al. 1998) and magnetic energy storage in active regions (Wang et al. 2009). 



2.5. Combining Logistic Growth and Diffusion Models 

In order to have compatible models for the exponential growth phase, we adopt the same exponential 
growth phase of the logistic model (Eq. 15) for the diffusion models, so that the spatio-temporal evolution 
r{t) differs only in the deceleration phase for the various models shown in Fig. 1. We choose the transition 
time between the acceleration and deceleration phase at t = ti in the formulation of the logistic growth curve 
(Eq. 15), which has the value ri and time derivative t'l, 

n = r{t = t,) = 2-V3^^ , (16) 

dr 2-4/3 Too ri 
vi = -7r{t = h) = — = — . (17) 

at S TG OTG 

Defining the diffusion radius r{t) in the deceleration phase with a diffusion constant /^, 

r{t) = K{t-tkf/\ (18) 
which has the following value ri and time derivative vi at the transition time t = ti , 

n = r{t = h) = K{h - tkf^^ , (19) 

which have to match the boundary conditions ri and vi (Eqs. 16-17) of the exponential growth phase. From 
Eqs. (16)-(20) we obtain the constants ti^VoQ^TG^ 

^ri\2//3 



ti = tk 



Too = ri 2^/3 , (22) 

TG = ^ . (23) 

Thus the combined time profile r{t) with a smooth transition from the exponential growth (acceleration) 
phase to the diffusion (deceleration) phase is. 



r{t)= { 'ooLl + exp(-^)j ... .^.x ^24) 

n{t-tkY/'^ ioTt>ti 

which has four free parameters for fitting, where we treat (/^, t/., /3, ri) as independent paramters, while 
(^i5 '^cx)5 '^g) are constrained by the boundary conditions (Eqs. 21-23). For classical diffusion (/3 = 1) the 
number of free parameters reduces to three. 



-1/3 



for t < ti 
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3. OBSERVATIONS AND DATA ANALYSIS 

3.1. GOES and AIA Observations 

We select all solar flare events detected with the Geostationary Operational Environmental Satellites 
(GOES) and the Atmospheric Imaging Assembly (AIA) on the Solar Dynamics Observatory (SDO) (Lemen 
et al. 2012) above a threshold of the Ml.O class level (which includes M- and X-class events) during the 
first two years of the SDO mission. The selected time era starts when the first science data from AIA 
became available, 13 May 2010, and ends on 31 March 2011 when we started the data analysis. The selected 
time period contains a total of 155 solar flare events larger than Ml.O class, including 12 events larger than 
XI. class. The flare time intervals, listed in Table 1, were taken from the official GOES fiare catalog. We 
extracted AIA images in all 6 coronal wavelengths during the GOES ffare time intervals, with AIA data 
being available in 100% of the cases, but restricted our analysis to the 335 A filter (Fe XVI), which appears 
to be the most suitable filter to probe the high-temperature component of these largest flares, since the other 
filters that are sensitive to high temperatures (94 A and 131 A) have a double response to cooler plasma of 
T < 1.0 MK that is not yet well- calibrated (due to the incompleteness of the CHIANTI atomic data base 
for Fe lines at temperatures < 10^ K). All AIA images have a cadence of At = 12 s and a pixel size of 
Ax = ^ 435 km, which corresponds to a spatial resolution of 2. 5 Ax = 1.5'' 1100 km. We normalized 
ah AIA 335 A images by the exposure time. The total number of analyzed images amounts to 11,767 AIA 
335 A images, which averages to 76 images per flare event and an average flare duration of 910 s. 



3.2. Data Analysis Method 

In order to quantify the spatio-temporal evolution of flares we measure the flare area above some 
threshold flux level, which was chosen in the AIA 335 A filter at a constant level of Fthresh = 100 and 200 
DN/s. This is well below the median fiux value of Fmed{t = 0) = 518 DN/s and mean of < F(t = 0) >= 
687 =b 645 DN/s of the maximum flux values of the EUV images at the flare start time, deflned by GOES. 
Thus the instantaneous flare area a(t) at time t is deflned by the number of pixels that have a flux F in 
excess of the threshold Fthresh , 

a{ti) = N[F^^y{ti) > Fthresh] , (25) 

while the time-integrated flare area A(t) is the combined area of all spatially overlapping instantaneous flare 
areas a(< t), 

A(ti) = a{t = 0) e a{t = 1) e ... e a{t = ti) , (26) 

where the symbol indicates a logical OR-function between the pixels contained in each instantaneous flare 
area before time ti. Thus, the time-integrated flare area A(t) is a monotonically increasing quantity that 
contains all flux pixels that exceeded the threshold Fthresh at any time before a given time t^, i.e., during 
the flare time interval < t < t^. 

The spatial scale r(t) of each (time-integrated) flare area A{t) is then deflned by the radius of an 
equivalent circular area, 

T{t) = ./mh ■ (27) 

The advantage to use the time-integrated flare area A(t) over the instantaneous flare area a(t) is the ro- 
bustness against temperature effects such as conductive or radiative cooling, which can decrease the flux of 
flaring pixels below the flux threshold after some time. However, we tested our method for both options and 
found compatible results in the cases of flares with no rapid cooling. 
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3.3. Examples of Analyzed Flares 

In Fig. 2 we show an example of the data analysis of a single event (#28), out of the 155 analyzed 
events. The GOES time profiles are shown on a linear flux scale (Fig. 2, top panel), with a start time of 
2011 March 07, 19:43 UT, peak time 20:12 UT, end time 20:58 UT, and GOES class M3.7 (= 3.7 x 10"^ W 
m~^), according to the official NOAA flare catalog. The flare duration is over an hour (T = 4500 s). 

At this point it might be appropriate to recall the definition of the flare duration used in the NOAA 
flare detection algorithm: The event starts when 4 consecutive 1-minute X-ray values have met all three of 
the following conditions: (i) All 4 values are above the Bl threshold; (ii) All 4 values are strictly increasing; 
(iii) The last value is greater than 1.4 times the value that occurred 3 minutes earlier. The peak time is 
when the flux value reaches the next local maximum. The event ends when the current flux reading returns 
to half of the peak value (|http: / / www.ngdc.noaa.gov/stp / solar/ solarfiares.html) . 

We analyzed the same flare time interval from AIA 335 A data, which amounts, with a cadence of 
AT = 12 s, to a total of Nt = T/AT = 4500/12 = 375 images. The total solar flux Fsun{t) during this 
time interval (Fig. 2, second panel) is shown along with the flux Ffiare{t) integrated over the flare area 
(approximately corresponding to the field-of-view shown in the bottom rows of Fig. 2). The agreement 
between the two flux profiles confirms the correct localization of the flare position, which is centered at the 
heliographic position N30/W48 for this event. 

The flare area A{t) is evaluated by counting the pixels above the threshold Fthresh = 100 DN/s, (or 200 
DN/s, respectivtly). The pixels of the time-integrated area combines all overlapping flare areas from earlier 
times t > to, since flare start. The radius r(t) = ^^ A(t)/7: of the equivalent circular (monotonically growing) 
fiare area A{t) is shown in the third panel of Fig. 2 (histogrammed curve), fitted with the model r{t) defined 
in Eq. (24) (shown with a solid curve in Fig. 2). The goodness-of-fit is simply evaluated from the average 
deviation normalized by the maximum flare radius, Vmax = max[r{t)] = r{tend)^ 

Ar/r™„. = <^o..W-r^o.eKt)> ^^^^ 

'^max 

which is Ar/rmax = 1-3% for the threshold of Fthresh = 100 DN/s, or Ar/rmax — 1-2% for the threshold 
of 

Fthresh — 200 DN/s, respectively. We determine the flare area for two different thresholds, in order to 
test the uncertainty and sensitivity of the best-fit model parameters on the threshold value. In this flare we 
obtain a diffusion coefficient of k, = 117 ± 10 pixel s~^/^ and a diffusion index of /3 = 0.47 ± 0.08, which 
demonstrates a relatively small uncertainty caused by the chosen two fixed flux thresholds (see values given 
in Fig. 2 and listed in Table 1). 

Five snapshots of the time- integrated flare images are shown in the bottom panels of Fig. 2, rendered as 
preflare-subtracted flux (Fig. 2, fourth row), as well as highpass-filtered flux (Fig. 2, bottom row), showing 
the flare loop fine structure with enhanced contrast. The contours of the flare area A{t) are also shown as 
black contour lines in the fiare snapshot images. The flux A{t = 0) is near zero at the beginning of the time 
series, indicating that the flare start time (defined by GOES) corresponds to an early phase when the 335 A 
flux was near the chosen threshold. 

The data analysis is also shown in condensed form for a selection of another 48 events (out of the 155 
analyzed events) in Fig. 3, where we show the spatio-temporal evolution of the GOES 1-8 A flux fcoEsif) 
(dotted curve), the EUV flux fAiA{t) of the 335 A filter (dashed curve), and the spatial evolution of the 
radius r{t) of the equivalent circular flare area (histograms in Fig. 3), along with the best fit of the diffusion 
model (Eq. 24) for the two flux thresholds of Fthresh = 100 and 200 DN/s (solid curves in Fig. 3), as well 
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as a (prefiare-subtracted) time-integrated AIA 335 image at the end of the flare (right panels in Fig. 3). 
The examples show the enormous variety of morphological shapes of flare areas, which generally consist of 
groups or arcades of post flare loops (in different projections) with highly fractal substructures. We show 
also the contours of the flnal time-integrated flare area above a threshold of Fthresh = 100 DN/s in Fig. 3 
(right panels), which generally encompasses a contiguous flare area. The time proflles reveal that the GOES 
flux generally peaks earlier than the EUV flux, although they start at similar times. The evolution of the 
flare radius r{t) shows often a flrst expansion triggered by a precursor flare. In such cases, our flts apply to 
the main phase, while the area associated with the precursor is subtracted. Note the variety of /3-values in 
different flares, ranging from the logistic limit /3 > 1 to the classical diffusion limit of /3 ~ 1. The fits of the 
diffusion model r{t) shown in Fig. 3 illustrate the adequacy of the functional form (Eq. 24) that fits the data 
always within a few percent accuracy, as well as the consistency of results obtained from different thresholds. 



3.4. Statistical Results 

The statistics of the observed parameters is visualized in form of histogram distribution functions shown 
in Fig. 4, as well as in form of correlations between pairs of parameters in Fig. 5. The best-fit parameters L, 
/3, and qfu = Ar/vmax are also listed in Table 1 for each flare, and the ranges, means, standard deviations, 
and medians are summarized in Table 2. 

The length scale L, deflned as the spatial radius of the flare area after subtraction of possible preflare 
areas, i.e., L = \/r{tendY ~ ^{^startY^ was found in a range of L = 5 — 50 Mm. The flare time scale T, 
which we deflne from the GOES flare start time to the peak time, e.g., T = tpea/c = t starts has a range of 
T = 120 — 8460 s (i.e., from 2 min to over 2 hours). The GOES flux ranges for this sample of M and X-class 
flares from Fgoes = 10"^ to 6.9 x 10"^ W m-^. The AIA 335 A flux is found in the range of Faia = 7 x 10^ 
to 2.5 X 10^ DN/s. Then we measured also the average expansion velocity v = Vmax/ {t start — tend) = 3.3 — 103 
km s~^, as well as the maximum velocity during the flare time interval, Vmax — 

8 - 550 km s-^ These 

parameters {L^T, FgoeS: FAiA^v^Vmax) have approximate powerlaw distributions (Fig. 4, left and middle 
column) . 

The best-flt model parameters have nearly Gaussian distributions (Fig. 4, right column), which can 
be characterized by the mean and standard deviations: the diffusion coefficient = 56 ± 24 km s~^/^, the 
diffusion index f3 = 0.53 ±0.27, and the goodness-of-flt Ar/r^ = 2.2 ±0.7%. Thus, the diffusion index, which 
is found in a range of /3 = 0.04 — 1.35, covers the entire range from logistic growth (/3 ~ 0), to sub-diffusion 
{f3 ^ 1/2), and classical diffusion {(3^1), but extends much rarer into the regime of super-diffusion or Levy 
flights {/3 ~ 3/2). No case was found that approaches the limit of linear expansion (/3 = 2). 



3.5. Parameter Correlations 

Parameter correlations can reveal physical scaling laws. In Fig. 5 we show scatterplots between pairs of 
observed parameters (L, T, Fgoes ^ Faia, ^, ^max, i^, Let us flrst discuss how these parameters depend 
on the fundamental length scale L of flares (Fig. 5). The flare duration T shows only a loose correlation 
with the length scale L, i.e., T oc i;,2.o±i.2 (^pg^j^el a), or L oc T'O-Sio.s (^pg^j^gi Weak correlations are also 
found for the GOES flux Fgoes (F) (panel b), the mean expansion velocity v{L) (panel c), as a function of 
the length scale L. The strongest correlations are found between the AIA 335 A flux and the length scale, 
i.e., Faia{L) oc l2.4±o.5 j^p^nel e), or with the volume Faia{V) oc yO-8±o.2 (jy^^^i g)^ ^nd the diffusion 
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coefficient tv{L) ex i^O-9±o i (panel f). The former good correlation can be explained if the EUV emission is 
proportional to the total flare volume, in which case we expect Faia oc V oc L^. The latter good correlation 
can be explained for sub-diffusion and logistic growth, where the diffusion equation L oc K>{t — ti)^/^ shows 
only a weak dependence on the time duration due to the small value of the power index /3/2 <C 1, which 
yields almost a proportionality of L oc a^. 

Correlating all other parameter combinations we found five more cases that are worthwhile to mention, 
shown in the lower half of Fig. 5. The mean expansion velocity scales almost reciprocally with the flare 
duration, Vmax{T) oc T^-i oiO-^ (panel h), which is consistent with a maximum speed that is independent of 
the length scale, so that the mean length scale < L > is a constant (with uncorrelated scatter) and yields 
"^max =< L > /T (X . The maximum velocity shows also a good correlation with the mean velocity, 
i.e., Vmaxi"^) OC ^^-^^O-^ (panel k). There is also an expected trend of correlated EUV and soft X-ray fluxes, 
i.e., Faia oc ^^(^0^5^ (panel i). What is most interesting that comes out of this study, is that the diffusion 
coefficient is strongly correlated with the EUV flux, i.e., k oc F^j^^'^ (panel 1), as well as with the soft X-ray 
flux, i.e., K oc F^Q%Y (panel j). 

4. DISCUSSION 

In this study we measured spatial (L) and temporal scales (T) in solar flares and investigated for the 
flrst time their spatio-temporal evolution L{t). We quantifled the spatio-temporal evolution with a general 
diffusion equation that is quantified in terms of a diffusion coefficient and a diffusion powerlaw index /3 
(Eq. 24). The theoretical range of the diffusion powerlaw index, /3 = — 2 includes the limit of logistic growth 
(/3 0), sub-diffusion (/3 ^ 1/2), classical diffusion (^ = 1), and super-diffusion or Levy fiights {/3 = 3/2). 
In the following we discuss the consequences of the results in the context of self-organized criticality models 
(section 4.1), diffusion processes in the photosphere (section 4.2), in the corona (section (4.3), and in solar 
flares (section 4.4). 

4.1. Self- Organized Criticality Models 

The size distribution of solar flare parameters (energy, peak count rates, durations) follow all powerlaw 
distribution functions that have been interpreted in terms of a slowly-driven nonlinear system in the state 
of self-organized criticality (SOC) (Lu and Hamilton 1993). A quantitative model that predicts the values 
of the powerlaw slopes and the associated correlations of physical parameters has been recently proposed, 
was tested with numerical simulations of cellular automaton models, and compared with solar flare statistics 
(Aschwanden 2012). The most fundamental parameters of SOC models are the spatial size (e.g., length 
scale L, area A, volume V) and the time duration T. In SOC models applied to astrophysical observations, 
generally a proportionality is assumed between the spatial volume V and the radiated fluence which can 
be justified by the linear relationship between the column depth / dz and emission measure EM oc J n^dz 
for optically-thin EUV and soft X-ray emission. Consequently, the instantaneous fractal volume of a SOC 
avalanche {dVs/dt) corresponds to the energy dissipation rate {dE/dt) or flux F, and the peak energy 
dissipation rate {dE / dt)max to the peak flux P. 

While classical diffusion (L oc T^/^) was adopted in the statistical fractal-diffusive avalanche model 
of a slowly-driven self-organized criticality system (FD-SOC) (Aschwanden 2012), based on empirical evi- 
dence from cellular automaton models, we generalize the model for anomalous diffusion and for wavelength- 
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dependent flux- volume scaling law Fx cx {dVs/dt)^ in Appendix A, which predicts slightly different correla- 
tion coefficients and powerlaw slopes of the distribution functions. In 3D-space (5 = 3), for which a mean 
fractal dimension of D3 = 2.0 is predicted, and using the measured range of the diffusion index P ~ 0.1 — 1.0 
(Fig. 4) and 7 = 0.8 (Fig. 5, panel g), the following powerlaw slopes are predicted by the generalized FD-SOC 
model (Appendix Eq. A14), 

aL =3.0 

aT = 1 + (S' - l)/3/2 ^ 1.1 - 2.0 

ap = 1 + (5' - l)/Dsj = 2.25 . (29) 

ap =l-h(5-l)/5'7=l-83 

ap =l^{S-l)/{Dsj^ 2/ 13)^1.1 -1.6 

For the size distributions of flare length scales L we expect in the framework of the FD-SOC model ap = 
while our observations show ap = 2.0 (Fig. 4, first row, left panel). For the size distribution of fiare durations 
T, classical diffusion predicts a slope of ap = 2.0 and sub-diffusion (with (3 = 1/2) predicts ap = 1.5, while 
our observations show ap ^ 2.2 (Fig. 4, second row, left panel). For the GOES flux Fq and AIA 335 flux 
Fa we measure powerlaw slopes of = 1.92 and ap = 1.34 (Fig. 4), while the FD-SOC model predicts 
ap = 2.0. In a much larger sample of over 300,000 GOES flares, a value of ap = 1.98 ±0.11 was found 
(Aschwanden and Freeland 2012), which is fully consistent with the FD-SOC model for classical diffusion 
{(3 = 1). Considering the small statistical sample of 155 flares analyzed here, the obtained powerlaw slope 
values seem to be not inconsistent with the FD-SOC model, given the relatively small range of less than a 
decade over which a powerlaw slope could be fltted. 

Alternatively, we can compare the observed parameter correlations (Fig. 5) and compare with the 
predictions of the SD-SOC model (Appendix A). The most fundamental relationship for anomalous diffusion 
predicts L oc T^/^, whith amounts to L oc T^^^ for sub-diffusion {(3 = 1/2), while our linear regression flt 
yields L oc T'O-Sio.s j^p^g panel d), which has a large uncertainy due to the weak dependence and large 
scatter of L and T values. A tighter correlation is found between the EUV flux and the flare volume, i.e., 
Faia cx: yO-8±o.2^ which corroborates the standard assumption in SOC models that the observed flux or 
emission is approximately proportional to the emitting volume. 

In summary, the small statistical sample of 155 flare events does not allow us to determine the statistical 
scaling laws (or correlations) between spatial L and temporal scales T or the powerlaw distributions N(L) 
and N{T) of spatial and temporal parameters with suflicient accuracy, and to test SOC models based on their 
statistical distributions. However, our study can quantify the spatio-temporal evolution L{t) for individual 
flare events very accurately, which could potentially be compared with the spatio-temporal evolution L{t) of 
cellular automaton simulations to test SOC models quantitatively. This could clarify whether SOC cellular 
automaton avalanches based on (thresholded) next-neighbor interactions can be statistically described by 
anomalous or classical diffusion processes. Diffusion-dominated dynamics (Bregman 2007, 2008) or diffusion 
entropy (Grigolini et al. 2002) has been applied to SOC avalanche models elsewhere. 



4.2. Diffusion Processes in the Photosphere 

Magnetic fluxtubes that emerge through the photosphere and build footpoints of coronal magnetic 
loops as well as postflare loops, are buffeted around by the random walk dynamics of subphotospheric 
magneto-convection. Measurements of the random motion of granules and supergranulation features have 
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been characterized with a diffusion equation, 

< >= ADt , (30) 

where the diffusion coefficient D is related to our definition V< > = A^t^s (Eq. 18) hjD = Thus, 
our measurements in solar flares with /i: = 56 =b 24 km s~^/^ would correspond to a range oi D ^ 300 — 1600 
km^ s~^. In comparison, photospheric random motion has been determined to D = 110 km s~^ in strong- 
field cores of active regions, and D = 250 km s~-^ in the surrounding area (Schrijver and Martin 1990), 
D = 200 - 400 km s"^ from the motion of magnetic elements (Mosher 1977), D = 120 - 230 km s"^ 
from Quiet Sun magnetic patterns (Komm et al. 1995), which appears substantially lower than in flares. 
However, higher values have been inferred from the magnetic polarity reversal during a solar cycle {D = 1100 
km s-^ Leighton 1964, 1969), from meridional flows {D = 600 km s-^ Wang et al. 1989; Zirin 1985), or 
from mottles {D = 600 km s~^ corrected for unresolved small-scale elements; Schrijver et al. 1996). In 
analogy to our fractal-diffusive SOC model (Aschwanden 2012), the dispersion of magnetic fleld elements 
across the solar surface was thought to be a diffusion process on a fractal geometry (Lawrence 1991). The 
diffusion of magnetic elements across the solar surface was actually found to have a 2-D fractal dimension 
of D2 = 1.56 ± 0.08 and an anomalous diffusion index of (3 = 0.25 ± 0.40, similar to our result of preferential 
sub-diffusion (/3 = 0.53 ± 0.27). In contrast, Ruzmaikin et al. (1996) found super-diffusion in photospheric 
random motion. A super-diffusive regime with /3 « 1.48 — 1.67 was also measured from the proper motion 
of bright points (Abramenko et al. 2011). Diffusion in the photosphere is also thought to be responsible for 
magnetic flux cancellation (Litvinenko 2011). 

In summary, most diffusion processes observed in magnetic features in the photosphere are generally 
slower {D ^ 100 — 600 km s~^) than what we observe in flares {D ^ 300 — 1600 km^ s~^), and thus are 
unlikely to explain the observed spatial diffusion during flares. 



4.3. Diffusion Processes in the Corona 

The magnetic fleld evolution in the solar corona is usually described by the induction equation in the 
framework of ideal magnet ohydro dynamics (MHD), 

1^ = V X (v X B) + r/V'B , (31) 

with T] = c^/ Ana being the magnetic diffusivity. The flrst term on the right-hand side is called the convective 
term, while the second term is called the diffusive term. Depending on the value of the Reynolds number 
Rm, which gives the ratio of the convective term (oc vqBo/Iq) to the diffusive term (oc tjBq/Iq), 

Rm = ^ , (32) 

the induction equation can be approximated in the two limits by 

r/B 

V X (v X B) for i?^ > 1 (33) 

^ « T^V^B for Rm <C 1 (34) 

The plasma in the solar corona is close to a perfectly conducting medium (with a high Reynolds number 
Rm ~ 10^ — 10^^), so that approximation (Eq. 33) with R^ ^ 1 applies, while the diffusive limit (Eq. 34) 
with Rm <^ 1 is not relevant. Thus, cross-fleld diffusion is generally inhibited in the solar corona. 
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While cross-field transport is generally inhibited in low plasma- /3 parameter regions, it can be enabled 
by gyro-orbit perturbations due to local magnetic field fluctuations, a process called Rechester-Rosenbluth 
diffusion (Galloway et al. 2006), but this process is estimated to operate on time scales of days, and thus 
is way too slow to explain the observed anomalous diffusion in flares. A scenario with small-scale random 
footpoint motion that injects energy into the corona by twisting and braiding was invoked as a coronal heating 
mechanism, where helicity conservation leads to hyper- diffusion (VanBallegooijen and Cranmer 2008), but 
this is also driven by photospheric diffusion speeds, and thus not fast enough to drive solar flares. 

4.4. Diffusion Processes in Solar Flares 

Let us discuss first diffusion processes in the prefiare phase. Resistive diffusion of force- free magnetic 
fields has been considered to develop infinite field gradients and to trigger the eruption of solar fiares (Low 
1973a,b; 1974b) and acceleration of particles (Low 1974a). On the other side, photospheric diffusion is 
thought to reduce the non-potential magnetic (free) energy that is built up by the photospheric shearing mo- 
tion (Wu et al. 1992). The eruption of a coronal mass ejection (CME) was simulated with MHD simulations, 
where slow turbulent diffusion of the footpoints of coronal magnetic field lines leads to a loss of equilibrium 
and drives the eruption (Amari et al. 2003). However, all these diffusion processes occur before a flare, and 
thus do not explain the rapid diffusion during fiares. 

What happens during fiares? The driving instability of solar fiares is generally linked to the dynamics 
of magnetic reconnection processes. Diffusion of particles is most efficient in the X-point of collisionless 
magnetic reconnection regions, where the magnetic field drops to zero. This diffusion region is subject 
of numerous theoretical and simulation studies (e.g., see review by Hesse et al. 2011). The duration of 
a magnetic reconnection process is approximately the Alfvenic transit time across a coronal reconnection 
region, which is of order Ti = L/'u^ « 10 s (for L = 10, 000 km and va = 1000 km s~^). While such a single 
magnetic reconnection episode (Sweet-Parker-type or Petschek-type) can only explain very short single-loop 
fiares, large fiares (with a range of durations T = 120 — 8460 s observed here) require a chain reaction 
of magnetic reconnection events, perhaps in the order of N = T/Ti = 12 — 850 reconnection episodes. 
Evidence for such multi-reconnect ion events has been demonstrated for numerous large flares, the most 
prominent case being the Bastille-Day flare, where over N > 200 postfiare loops have been traced, probably 
each one being a remnant of a local magnetic reconnection process (e.g., Aschwanden and Alexander 2001), 
forming a multi-loop postfiare configuration (Hori et al. 1998). Theoretical models of multi-reconnection or 
intermittent, unsteady, bursty reconnection have also been formulated in terms of the tearing mode instability 
and coalescence instability, and combinations of both (e.g., Furth et al. 1963; Sturrock 1996; Kliem 1990, 
1995; Leboef et al. 1982; Tajima et al. 1982; Karpen et al. 1995; Kliem et al. 2000; Drake et al. 2009). 
The question is now how the time evolution of such an intermittent chain reaction of magnetic reconnection 
episodes ties into our simplified diffusion evolutionary model. 

Perhaps our measurements of the average expansion speed can give us a hint about the physical mecha- 
nism of the diffusion process. We measured average speeds in the range oiv = 3 — 103 km s~^, and maximum 
speeds of Vmax = 8 — 550 km s~-^. If we assume typical plasma temperatures of T 10 — 30 MK for the 
observed large (M and X-class) fiares, we expect sound speeds of Cs ~ I^QV^mk ^ 500 — 900 km s~^, 
which are about 1-2 orders of magnitude faster than the observed speeds. Therefore, the diffusive speed we 
observe in these large fiares cannot just be the propagation speed of upfiowing plasma into and along fiare 
loops initiated by the chromospheric evaporation process, which occurs approximately with sound speed. 
There are some additional delays that intervene between the subsequent filling of individual loops in a flare 
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arcade. A model where subsequent flare loops are triggered by a slow (acoustic) wave that propagates at 
some angle (25° — 28°) along the arcade was proposed by Nakariakov and Zimovets (2011) and tested by the 
timing of the hard X-ray footpoints (Inglis and Dennis 2012), but the predicted hard X-ray pulse periods 
did not agree with the observed ones. Nevertheless, the mean footpoint speed found in the analyzed three 
flares {v = 5 — 60 km s""*^; Inglis and Dennis 2012) is similar to our observed range {v = 3- 103 km s-^), 
which corroborates the congruence between the magnetic reconnection path and the hard X-ray footpoint 
path. Similar propagation speeds of the hard X-ray footpoints were also measured along an arcade of flare 
loops by Grigis and Benz (2005), i.e., a mean velocity of t' = 63 — 65 km s~-^ parallel to the arcade, with 
a peak speed of about v = 110 km s~^ for 2 minutes. Hence we can interpret our diffusion speed as mean 
propagation speed of subsequently triggered magnetic reconnection sites. 

The measurement of the spatio-temporal evolution during solar flares should also ameliorate estimates 
of the particle number problem (Brown and Melrose 1977), which requires the knowledge of the flare volume 
and its temporal change and spatial propagation path. 



4.5. Fractal Diffusion, Sub-DifFusion, or Logistic Growth ? 

We combined the different anomalous and classical diffusion models into a single framework (Eq. 24) 
that has only one single parameter (/3) that discriminates between the different models: /3 = is the limit 
of logistic growth, ^ = — 1 is the sub-diffusive regime, (3 = 1 is classical diffusion or Brownian motion, 
P = 1 — 2 is the super-diffusive regime or Levy flights, and /3 = 2 is the limit of linear expansion. Our 
results yielded a range of /3 = 0.04 — 1.35 that covers three of these regimes, but the majority is found in the 
sub-diffusive regime with a mean and standard deviation of /3 = 0.53 =b 0.27. What does this mean and why 
is it different from the classical diffusion model based on random walk? 

Brownian motion or classical diffusion is based on homogeneous isotropic expansion. In solar flares, as 
well as in cellular automaton models, however, the medium is inhomogeneous (which can be characterized 
by a fractal dimension). The solar corona, moreover, is highly anisotropic due to the structuring by the 
magnetic fleld (in the low plasma- /3 regime), which also applies to flares, except to a reduced degree in the 
diffusion regions of magnetic X-point reconnection conflgurations. Let us assume the extreme case where 
plasma can only flow in one direction along magnetic fluxtubes. The anisotropic volume increase will scale 
as V{t) oc r(t), instead of V{t) oc in an isotropic medium. The observable area A{t) will then grow only 
in one direction (say in x-direction) with a constant width (say A^/), 

Aiso(t) = Xiso{t)Ay , (35) 

and the derived flare expansion radius r{t) = y^A{t)/7T will scale as 

r{t) = ^J^^ oc xUt)'^' oc t^l' , (36) 

for a classical diffusion process in one dimension, i.e., Xiso{t) oc t^/^. Thus by deflning anomalous diffusion 
by r{t) oc t^/^ we end up with a diffusion index of /3 = 1/2, which is indicative of the sub-diffusive regime. 
Now, although the electron and ion diffusion region in magnetic reconnection X-points allow 2-D or 3-D 
diffusion, our observational result of sub-diffusion with P ^ 1/2 suggests that the triggering of subsequent 
magnetic reconnection episodes occurs along anisotropic 1-D paths and dominates the overall diffusion of 
the flare process, regardless of the near-isotropic diffusion regions in X-points. 
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What does the extreme limit of logistic growth mean, where the diffusion comes to a halt at some 
time near the peak time of the flare and does not expand further. Apparently, the energy release process 
approaches regions with strong magnetic fields that are stable and not prone to magnetic reconnection, so 
that the chain reaction of reconnection episodes does not propagate further, while energy release in the 
previous area still continues. It is like a domino chain reaction that hits the wall of a room. The domino 
effect can still propagate in the unstable regions inside the room, but cannot expand outside of the room. 
Alternatively, the progressing flare could also reach the boundary of an active region, where only weak-field 
regions of the Quiet Sun are available which harbour much less non-potential magnetic energy than the 
inside of active regions. 

In summary, the diffusion index {3 tells us some interesting information about the anisotropy and bound- 
aries of magnetic topologies in solar flare regions. 

5. CONCLUSIONS 

We analyzed the spatio-temporal evolution in the 155 largest solar flares (M and X-class) observed by 
GOES and AIA/SDO during the first two years of the SDO mission. We fitted the radial expansion r{t) of 
flare areas detected in the 335 A filter above some threshold level with a generalized diffusion model that 
includes the classical diffusion, anomalous diffusion, and the logistic growth limit. The major results and 
conclusions are: 

1. The flare area in all events can be fitted with a radial expansion model r{t) that consists of an initial 
acceleration phase with exponential growth and a deceleration phase that follows anomalous diffusion, 
r{t) oc A^t^/^, with /3 = 0.53 ±0.27 mostly falling into the sub-diffusive regime. The most extreme cases 
range from logistic growth (/3 = 0.04) to super-diffusion (/3 = 1.35). The sub-diffusive characteristics is 
likely to reflect the anisotropic propagation of energy release in a magnetically dominated plasma. The 
limit of logistic growth indicates the times when the boundaries of energy release regions are reached. 

2. The diffusion coefficient = 53 dz 23 km s~^/^, which corresponds to an area diffusion constant of 
D ^ 300 — 1600 km^ s~^ is found to be significantly faster than diffusion processes measured in the 
photosphere {D ^ 100 — 600 km s~^) and cannot be explained with crossfield diffusion in the corona, 
which is strongly inhibited by the low plasma- /3 parameter. 

3. The average diffusion speed during flares is measured in the range of v ^ 5 — 100 km s~^, with 
maximum speeds of Vmax ~ 10 — 500 km s~^, which is slower than the sound speed of Cs ~ 500 — 900 
km s~^ expected in flares with temperatures of Te ~ 10 — 30 MK, but is compatible with the hard 
X-ray footpoint motion along the neutral line during flares, and thus is likely to represent the mean 
propagation speed of subsequently triggered magnetic reconnection sites. 

4. The fractal- diffusive self-organized criticality model (FD-SOC) describes the diffusive flare progression 
in a fractal geometry and predicts powerlaw distributions with slopes of aL = 3.0 for length scales L, 
aT = 1.1 — 2.0 for time scales T, and ap = 1.83 for peak fluxes P, based on the observed diffusion 
index range of (3 ^ 0.1 — 1.0, the observed flux- volume scaling of F335 oc (dV/dt)^-^, in a 3D-space 
geometry with fractal dimension Ds = 2.0. Our observations yield q^l = 2.0, ar = 2.2, and ap = 1.34, 
which is not inconsistent with the FD-SOC model, given the large uncertainties of the powerlaw slopes 
determined in a small statistical sample. Nevertheless, the evolutionary fits r(t) oc t^/^ of individual 
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flares confirm the diffusive property of tiie FD-SOC model for tiie case of solar flares with a high degree 
of accuracy, which has been previously predicted based on cellular automaton simulations. 

Thus, the major accomplishment of this study is the experimental proof of a diffusive flare expansion 
process that has been predicted by the fractal-diffusive self-organized criticality (FD-SOC) model. How- 
ever, while the classical FD-SOC model assumed classical diffusion (/3 = 1), the observational results from 
AIA/SDO revealed mostly sub-diffusion (/3 < 1), down to the limit of logistic growth (/3 ^ 0). Future studies 
may extend the statistics of spatio-temporal parameters and allow us to test the predicted size distributions 
of SOC models with higher accuracy and establish physical scaling laws for flares that are based on the 
fundamental geometric parameters of space and time. 
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Appendix A: Generalized Fractal-Diffusive Self- Organized Criticality Model 

We generalize the fractal- diffusive self-organized criticality model (FD-SOC) described in Aschwanden 
(2012) in two respects, (1) regarding anomalous diffusion, and (2) for arbitrary (wavelength-dependent) flux- 
volume scaling. In the original FD-SOC model, classical diffusion was assumed, i.e., a scaling of L oc T^/^ 
between the length scale L and time duration T, which we generalize for anomalous diffusion (with diffusion 
index /3), 

L oc . (Al) 

which includes also the case of classical diffusion {(3 = 1). The second generalization is the the scaling of an 
observed flux Fx in a given wavelength A to the emissive volume rate of change dVs/dt^ 

where 7 = 1 was assumed to be unity in the original FD-SOC model, but can now be wavelength-dependent 
and accomodate some astrophysical scaling law of the underlying emission and absorption process at wave- 
length A. 

The FD-SOC model describes the inhomogeneity of SOC avalanches with a fractal structure with Haus- 
dorff dimension Ds in Euclidean space with dimension S. The instantaneous avalanche volume dVs/dt is 
then characterized by 

and the mean fractal dimension Ds can be estimated from the arithmetic mean of the minimum and maxi- 
mum possible value, 

yielding a mean value of D3 = 2.0 for the 3D Euclidean space (S* = 3). Combining Eqs. (A1-A3) yields then 
the correlation of the instantaneous energy dissipation volume (dVs/dt) and the time duration T, 

oc (^) ^ oc L^^^ oc T^s7/3/2 . (^5) 
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This flux quantity Fx is generally strongly fluctuating, and the peak dissipation rate P\ = max{F\) can be 
estimated from the maximum fractal dimension Ds $ 



Px cc 



V dt 



{A6) 



A more steady quantity is the time-integrated volume of a SOC avalanche, which can be obtained from 
the time- integration of the instantaneous volume rate of change {dVs/dt), and be associated with the fluence 
Ex (or time-integrated flux E\ = j F\{t)dt)^ 



Ex oc r Fx{t) dt oc r ^ 
Jo Jo 



{A7) 



After we obtained the scaling between the parameters L,T, Fx, Px, and Ex, we can calculate the powerlaw 
slopes of their size distributions. The basic reference is the distribution function N{L) of spatial scales, 
which be reciprocal to the size in a slowly-driven homogeneous SOC system that is scale- free and where 
avalanches can occur in any size. The relative probability is then (Aschwanden 2012), 



N{L) oc oc L- 



All other frequency distributions can then be derived by substituting the corresponding parameter correla- 
tions given in Eq. (A1-A7), First we can calculate the occurrence frequency distribution of avalanche time 
scales, by using the diffusive boundary propagation relationship L{T) oc T^/^ (Eq. Al), by substituting the 
variable T for L in the distribution N{L) (Eq. A8), 



N{T)dT = N{L[T]) 



dL 



dT 



dT oc T-^^HS-i)^/2] ^j. ^ 



Subsequently we can derive the occurrence frequency distribution function N{F) for the instantaneous energy 
dissipation rate F using the relationship Fx{T) oc T^sj0/2 ^5), 



N{Fx)dFx = N{T[Fx]) 



dT 



dFx 



(^10) 



the occurrence frequency distribution function of the peak energy dissipation rate P using the relationship 
relationship Px{T) oc T^"'/^/'^ (Eq. A6), 



N{Px)dPx = N{T[Px]) 



dT 



dPx 



dP,ocP-[^+(^-i)/^^ldP, 



(All) 



and the occurrence frequency distribution function N{E\) for the total energy E\ using the relationship 
ExiT) oc r(i+^s7/3/2) (Eq. A7), 



NiEx)dEx = NiT[Ex]) 



dT 



dEx 



dEx oc E7^'+^'-'^/^'''^'^+'/^'>Ue . 



(AU) 



This derivation yields naturally powerlaw functions for all parameters L, T, Fx, Px, and Ex, which are 
the hallmarks of SOC systems. In summary, if we denote the occurrence frequency distributions N{x) of a 
parameter x with a powerlaw distribution with power index ax , 



N{x)dx (XX dx , 



(AU) 
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we have the following powerlaw coefficients olx for the parameters x = L^T^F\^P\^ and 

aL = S 

ar = 1 + (5 - l)/3/2 

ap = 1 + (5 - l)/(i^57) • (^14) 

ap = 1 + (5 - 

ap =l^{S-l)/{Dsl^2/P) 

Note, that the powerlaw slopes q^l, Q^f and cep do not depend on the diffusion index /3, and thus are identical 
for classical or anomalous diffusion, while ap and ap depend on the diffusion index /3. 
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Table 1. Catalog of analyzed M and X-class flare events and best-fit model parameters: length scale 
L(Mm), diffusion coefficient k (km s~^/^), diffusion index /3, and goodness-of-fit qfu. 



iNr 


Observation 


Start 


Peak 


rL/na 


Duration 




IN OA A 


Heliogr. 


Length 


JJllT. 


JJlIT. 


r It 




date 


time 


time 


time 


l(s) 


class 


A TD 

ArC 


position 


L[NLin) 


coeff. K, 


index (3 




1 


2010-06-12 


00:30 


00:57 


01:02 


1920 


M2.0 


11081 


N23W47 


12 


65±25 


0.60±0.27 


2.9% 


2 


2010-06-13 


05:30 


05:39 


05:44 


840 


Ml.O 


11079 


S24W82 


8 


35ib 4 


0.81=b0.11 


1.8% 


3 


2010-08-07 


17:55 


18:24 


18:47 


3120 


Ml.O 


11093 


N14E37 


40 


107ib 4 


0.36±0.02 


2.2% 


4 


2010-10-16 


19:07 


19:12 


19:15 


480 


M2.9 


11112 


S20W26 


17 


86±18 


0.55±0.03 


2.6% 


5 


2010-11-04 


23:30 


23:58 


00:12 


2520 


Ml. 6 


11121 


S20E85 


12 


34± 


0.23±0.05 


1.9% 


6 


2010-11-05 


12:43 


13:29 


14:06 


4980 


Ml.O 


11121 


S20E75 


14 


38± 2 


0.58±0.06 


1.6% 


7 


2010-11-06 


15:27 


15:36 


15:44 


1020 


M5.4 


11121 


S20E58 


14 


68di 6 


0.53di0.01 


1.7% 


8 


2011-01-28 


00:44 


01:03 


01:10 


1560 


Ml. 3 


11149 


N16W88 


11 


30di 2 


0.12dz0.06 


3.3% 


9 


2011-02-09 


01:23 


01:31 


01:35 


720 


Ml. 9 


11153 


N16W70 


14 


46ibl2 


0.23ib0.09 


3.4% 


10 


2011-02-13 


17:28 


17:38 


17:47 


1140 


M6.6 


11158 


S20E05 


21 


74±12 


0.46±0.04 


2.6% 


11 


2011-02-14 


17:20 


17:26 


17:32 


720 


M2.2 


11158 


S20E04 


23 


66±11 


0.15±0.01 


2.2% 


12 


2011-02-15 


01:44 


01:56 


02:06 


1320 


X2.2 


11158 


S20W10 


33 


114± 3 


0.66±0.28 


2.9% 


13 


2011-02-16 


01:32 


01:39 


01:46 


840 


Ml.O 


11158 


S20W24 


12 


46± 5 


0.56±0.04 


2.5% 


14 


2011-02-16 


07:35 


07:44 


07:55 


1200 


Ml.l 


11161 


S19W29 


10 


41± 4 


0.66ib0.01 


2.5% 


15 


2011-02-16 


14:19 


14:25 


14:29 


600 


Ml. 6 


11158 


S21W30 


16 


41dz 4 


O.lOdiO.Ol 


1.1% 


16 


2011-02-18 


09:55 


10:11 


10:15 


1200 


M6.6 


11158 


S21W55 


13 


81ibl3 


0.57ib0.05 


4.2% 


17 


2011-02-18 


10:23 


10:26 


10:37 


840 


Ml.O 


11162 


N22E10 


14 


41± 4 


0.25±0.03 


2.6% 


18 


2011-02-18 


12:59 


13:03 


13:06 


420 


Ml. 4 


11158 


S20W70 


7 


47± 5 


0.81±0.03 


1.8% 


19 


2011-02-18 


14:00 


14:08 


14:15 


900 


Ml.O 


11162 


N22E10 


14 


40± 5 


0.17±0.00 


2.0% 


20 


2011-02-18 


20:56 


21:04 


21:14 


1080 


Ml. 3 


11162 


N22E10 


9 


48± 2 


0.81±0.09 


1.8% 


21 


2011-02-24 


07:23 


07:35 


07:42 


1140 


M3.5 


11163 


N14E87 


17 


77±17 


0.88dz0.12 


2.5% 


22 


2011-02-28 


12:38 


12:52 


13:03 


1500 


Ml.l 


11164 


N28E39 


20 


55ib 6 


0.28di0.02 


2.3% 


23 


2011-03-07 


05:00 


05:13 


05:19 


1140 


Ml. 2 


11164 


N23E47 


21 


71ibll 


0.58ib0.03 


2.2% 


24 


2011-03-07 


07:49 


07:54 


07:56 


420 


Ml. 5 


11165 


S18W75 


9 


38±15 


0.34±0.19 


2.7% 


25 


2011-03-07 


07:59 


08:07 


08:15 


960 


Ml. 4 


11164 


N27W46 


23 


63± 6 


0.18±0.02 


2.0% 


26 


2011-03-07 


09:14 


09:20 


09:28 


840 


Ml. 8 


11164 


S17W77 


5 


31±15 


0.92±0.33 


3.0% 


27 


2011-03-07 


13:45 


14:30 


14:56 


4260 


Ml. 9 


11166 


N11E21 


30 


72± 6 


0.42±0.04 


2.1% 


28 


2011-03-07 


19:43 


20:12 


20:58 


4500 


M3.7 


11164 


N30W48 


51 


117ibl0 


0.47ib0.08 


1.2% 


29 


2011-03-07 


21:45 


21:50 


21:55 


600 


Ml. 5 


11165 


S17W82 


10 


25± 1 


O.lOiO.Ol 


2.0% 


30 


2011-03-08 


02:24 


02:29 


02:32 


480 


Ml. 3 


11165 


S18W80 


9 


27± 3 


0.26±0.03 


2.0% 


31 


2011-03-08 


03:37 


03:58 


04:20 


2580 


Ml. 5 


11171 


S21E72 


32 


90± 5 


0.79±0.07 


1.7% 


32 


2011-03-08 


10:35 


10:44 


10:55 


1200 


M5.3 


11165 


S17W88 


19 


59± 9 


38±0 03 

\J • kJkJ —i—\J • KJkJ 


3.3% 


33 


2011-03-08 


18:08 


18:28 


18:41 


1980 


M4.4 


11165 


S17W88 


13 


37di 7 


0.52±0.15 


2.7% 


34 


2011-03-08 


19:46 


20:16 


21:19 


5580 


M1.4 


11165 


N23W43 


26 


60di 3 


0.50ib0.09 


2.0% 


35 


2011-03-09 


10:35 


11:07 


11:21 


2760 


Ml. 7 


11166 


N08W11 


21 


56ib 6 


0.63ib0.04 


3.2% 


36 


2011-03-09 


13:17 


14:02 


14:13 


3360 


Ml. 7 


11166 


S26W78 


12 


27± 7 


0.73±0.12 


2.8% 


37 


2011-03-09 


23:13 


23:23 


23:29 


960 


XI. 5 


11166 


N08W11 


34 


123±17 


0.56±0.03 


1.9% 


38 


2011-03-10 


22:34 


22:41 


22:49 


900 


Ml.l 


11166 


N08W25 


8 


42±10 


1.14±0.01 


2.3% 


39 


2011-03-12 


04:33 


04:43 


04:48 


900 


Ml. 3 


11166 


N07W35 


16 


53± 9 


0.35±0.03 


1.7% 


40 


2011-03-14 


19:30 


19:52 


19:54 


1440 


M4.2 


11169 


N16W49 


17 


54di 6 


0.26ib0.04 


3.4% 


41 


2011-03-15 


00:18 


00:22 


00:24 


360 


Ml.O 


11169 


N11W83 


9 


27ib 5 


0.17zb0.06 


2.7% 


42 


2011-03-23 


02:03 


02:17 


02:24 


1260 


M1.4 


11176 


S25W84 


18 


58ib 8 


0.45ib0.02 


1.4% 


43 


2011-03-24 


12:01 


12:07 


12:11 


600 


Ml.O 


11176 


S15E43 


15 


48± 5 


0.26±0.03 


2.3% 


44 


2011-03-25 


23:08 


23:22 


23:30 


1320 


Ml.O 


11176 


S12E26 


20 


59± 2 


0.39±0.06 


1.9% 


45 


2011-04-15 


17:02 


17:12 


17:28 


1560 


Ml. 3 


11190 


N13W24 


13 


43± 6 


0.42±0.04 


1.4% 


46 


2011-04-22 


04:35 


04:57 


05:14 


2340 


Ml. 8 


11195 


S17E42 


20 


60±12 


0.77±0.14 


2.0% 


47 


2011-04-22 


15:47 


15:53 


16:11 


1440 


Ml. 2 


11195 


S18E36 


20 


57± 5 


0.36di0.03 


1.5% 
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Table 1 — Continued 



Nr 


Observation 


Start 


Peak 


End 


Duration 


GOES 


NOAA 


Heliogr. 


Length 


Diff. 


DiflF. 


Fit 




date 


time 


time 


time 


T(s) 


class 


AR 


position 


L(Mm) 


coeff. K 


index /3 




48 


z011-05-z8 


21:09 


21:50 


22:01 


3120 


Ml.l 


11226 


COI T7"7n 

b21E70 


14 


37dz 5 


n TO-Ln no 

0.72±0.08 


1 no/ 
1.970 


49 


2011-05-29 


10:08 


10:33 


11:08 


3600 


Ml. 4 


11226 


b20E54 


24 


68± 6 


n 00 \ f\ 1 A 

0.93±0.14 


croy 
2.5% 


50 


2011-05-07 


06:16 


06:41 


06:59 


2580 


M2.5 


11226 


b22W53 


33 


85=b 7 


n on 1 n no 

0.30±0.02 


1 i::0/ 
1.5/0 


51 


2OII-O0-I4 


21:36 


21:47 


22:10 


2040 


Ml. 3 


11236 


ATI ATT^TT 

IN 14h;77 


20 


55zb 4 


n OT 1 n no 

0.37±0.02 


1 TO/ 

1.770 


52 


2011-07-27 


15:48 


16:07 


16:22 


2040 


Ml.l 


11260 


AT 1 OT7I 

JN19E38 


19 


63± 8 


0.55±0.15 


2.7% 


53 


2011-07-30 


02:04 


02:09 


02:12 


480 


M9.3 


11261 


AT 1 A XT' CT 

JN 14E35 


26 


1 on 1 on 

139±20 


n ^Jo_i_n ni 
0.52±0.01 


2.7% 


54 


2011-08-02 


05:19 


06:19 


06:48 


5340 


Ml. 4 


11261 


AT 1 <T\ X /■ 1 

JN17W12 


28 


74di 7 


n 1 n n cr 

0.58±0.05 


1 TO/ 

1.7% 


55 


2011-08-03 


03:08 


03:37 


03:51 


2580 


A /ri 1 
Ml.l 


11261 


ATI CATIT'OO 

i\ 15W28 


13 


37ib 4 


n /I _i_n n^j 
0.24±0.05 


aoy 
2.5/0 


56 


2011-08-03 


04:29 


04:32 


04:35 


360 


Ml. 7 


11263 


AT 1 r'TT^ 1 n 
IN 15ll;10 


21 


73zb 7 


n on 1 n ni 

0.29±0.01 


1 nO/ 
1.0%) 


57 


2011-08-03 


13:17 


13:48 


14:10 


3180 


M6.0 


11261 


ATI TTTS TOO 

JN17W30 


28 


78zb 3 


TO 1 /I 

0.72±0.04 


3.0% 


58 


2011-08-04 


03:41 


03:57 


04:04 


1380 


M9.3 


11261 


AT 1 r'W ro 

IN 15W38 


32 


92zbl6 


n 00 1 n no 

0.28±0.02 


1 nO/ 
1.9/0 


59 


2011-08-08 


18:00 


18:10 


18:18 


1080 


M3.5 


11263 


ATI CCTST'^JO 

IN low 52 


17 


57=b 6 


n /I ^7 1 n ni 

0.47±0.01 


1 ir:0/ 
1.5/0 


60 


2011-08-09 


03:19 


03:54 


04:08 


2940 


M2.5 


11263 


AT 1 >'7\ 1 T nc\ 

JN17W59 


12 


39di 5 


1 rr 1 n no 

1.35±0.03 


010/ 


61 


2011-08-09 


07:48 


08:05 


08:08 


1200 


X6.9 


11263 


ATI AWJClCy 
JN 14W59 


29 


95dil2 


n 1 n_i_n nn 
0.19±0.00 


3.2% 


62 


2011-09-04 


11:21 


11:45 


11:50 


1740 


M3.2 


11286 


AT 1 01 Tt ro A 

JN18W84 


10 


38ib 1 


n rrn i n i /i 

0.50±0.14 


3.3% 


63 


oni 1 no 
2011-09-05 


04:08 


04:28 


04:32 


1440 


Ml. 6 


11286 


ATI OTAT'OT 

IN 18W87 


13 


46dz 


n K K _Ln 1 
0.55±0.12 


nO/ 
2.0/0 


64 


2011-09-05 


07:27 


07:58 


08:06 


2340 


Ml. 2 


11286 


AT 1 OTS ro T 

IN 18W87 


7 


22d= 2 


1 00 1 n on 

1.33±0.29 


oO/ 

2.3%) 


65 


2011-09-Od 


01:35 


01:50 


02:05 


1800 


M5.3 


11283 


AT 1 OT 1 rOT 

JN 13W07 


21 


60zb 8 


T 1 r 

0.37±0.05 


0.9% 


66 


2011-09-05 


22:12 


22:20 


22:24 


720 


X2.1 


11283 


AT 1 /I "\ X n 

IN 14W18 


29 


87=b 6 


n 1 1 n no 

0.18±0.03 


/1 0/ 
3.4/0 


67 


2011-09-07 


22:32 


22:38 


22:44 


720 


XI. 8 


11283 


ATI AWJO'X 

JN 14W31 


36 


101dil5 


n 1 T-Ln no 
0.17±0.02 


1 00/ 
1.2/0 


68 


2011-09-08 


15:32 


15:46 


15:52 


1200 


M6.7 


11283 


AT 1 AWT A'\ 

IN 14W41 


27 


80dbl4 


n oo_i_n no 
0.32±0.08 


1 00/ 
1.8/0 


69 


2011-09-09 


06:01 


06:11 


06:17 


960 


M2.7 


11283 


AT 1 AWT AO 

JN 14W48 


16 


73dz21 


n ott_i_n on 

0.85±0.30 


1 00/ 
1.8/0 


70 


2011-09-09 


12:39 


12:49 


12:56 


1020 


Ml. 2 


11283 


ATI rTX TV r\ 

JN15W50 


15 


45zb 5 


0.35±0.01 


1 10/ 
1.1%) 


71 


oni 1 nn i n 
2011-09-10 


07:18 


07:40 


07:56 


2280 


Ml.l 


11283 


AT 1 /I T A JC2. A 

IN 14 W54 


15 


47d= 6 


n ni _Ln no 
0.91±0.02 


1 nO/ 
1.0%) 


72 


2011-09-21 


12:04 


12:23 


12:45 


2460 


Ml. 8 


11301 


AT 1 C TT* 

IN 15il;88 


6 


18=b 3 


n /1 1 n no 

0.48±0.03 


<:;0/ 
2.5/0 


73 


2011-09-22 


09:53 


10:00 


10:09 


960 


Ml.l 


11302 


ATo /IT X rrr rr 

JN24W55 


12 


52± 9 


0.87±0.04 


1 0/ 

3.1/0 


74 


oni 1 nn oo 

2011-09-22 


10:29 


11:01 


11:44 


4500 


XI. 4 


11302 


ATnoTZPon 

JN08E89 


41 


85ib 5 


n nn i n i o 

0.99±0.13 


1 00/ 
1.3/0 


75 


oni 1 nn oo 

2011-09-23 


01:47 


01:59 


02:10 


1380 


Ml. 6 


11295 


ATO /IT X J C A 

IN 24 W 54 


14 


59zbl0 


1 on 1 n no 

1.20±0.02 


oO/ 

2.3%) 


76 


oni 1 nn oo 
2011-09-23 


21:54 


22:15 


22:34 


2400 


Ml. 6 


11295 


ATI OTT* 

IN 12h;55 


17 


52d= 5 


n TO-Ln no 
0.72±0.08 


010/ 
2.1%) 


77 


oni 1 nn oo 

2011-09-23 


23:48 


23:56 


00:04 


960 
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1.8% 
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20:42 
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N13E52 


11 


47± 9 


0.55±0.03 


3.8% 


85 


2011-09-24 
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21:32 


540 


Ml. 2 
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N13E52 


5 


22±15 
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4.1% 


86 
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23:45 
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00:09 


1440 
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2.2% 


87 
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N22W87 


8 
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0.99ib0.00 


2.4% 


88 


2011-09-25 


04:31 


04:50 


05:05 


2040 


M7.4 


11302 


N13E50 


30 


90± 4 


0.58±0.08 


2.4% 


89 


2011-09-25 


08:46 


08:49 


08:52 


360 


M3.1 


11302 


N13E45 


11 


45± 6 


0.63±0.02 


1.7% 


90 


2011-09-25 


09:25 


09:35 


09:53 


1680 


Ml. 5 


11303 


S28W71 


10 


29± 


0.86±0.20 


4.1% 


91 


2011-09-25 


15:26 


15:33 


15:38 


720 


M3.7 


11302 


N13E44 


25 


72± 9 


0.19±0.02 


1.6% 


92 


2011-09-25 


16:51 


16:58 


17:09 


1080 


M2.2 


11303 


N12E41 


11 


39± 2 


0.51±0.06 


2.7% 


93 


2011-09-26 


05:06 


05:08 


05:13 


420 


M4.0 


11302 


N12E34 


17 


64±13 


0.43dz0.02 


2.9% 


94 


2011-09-26 


14:37 


14:46 


15:02 


1500 


M2.6 


11302 


N12E30 


28 


76± 5 


0.28ib0.00 


1.1% 
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iNi 


Observation 


otart 


r eaK 


Ji/na 


Duration 






Heliogr. 


Length. 


uin. 


uin. 


r It 




date 


time 


time 


time 


i {S) 


class 


A R 


position 


T /'A/TTY^^ 

ij^^ivim J 


coefF. Hi 


index /3 




95 


2011-09-28 


13:24 


13:28 


13:30 


360 


Ml. 2 


11302 


NllEOO 


16 


51ib 5 


0.25±0.02 


2.6% 


96 


2011-09-30 


18:55 


19:06 


19:15 


1200 


Ml.O 


11305 


N09E03 


18 


64ib 7 


0.67ib0.11 


2.6% 


97 


2011-10-01 


08:56 


09:59 


10:17 


4860 


Ml. 2 


11305 


N09W04 


24 


68± 6 


0.84±0.44 


2.2% 


98 


2011-10-02 


00:37 


00:50 


00:59 


1320 


M3.9 


11305 


N10W13 


18 


59± 8 


0.45±0.05 


1.6% 


99 


2011-10-02 


17:19 


17:23 


17:26 


420 


Ml. 3 


11302 


N10W55 


12 


41± 3 


0.27±0.00 


2.2% 


100 


2011-10-20 


03:10 


03:25 


03:44 


2040 


Ml. 6 


11318 


N18W88 


11 


28ib 6 


0.27zb0.14 


3.4% 


101 


2011-10-21 


12:53 


13:00 


13:08 


900 


Ml. 3 


11319 


N05W79 


12 


49di 


0.60±0.13 


2.5% 


102 


2011-10-22 


10:00 


11:10 


13:09 


11340 


Ml. 3 


11314 


N27W87 


50 


112ib23 


0.23ib0.02 


1.7% 


103 


2011-10-31 


14:55 


15:08 


15:27 


1920 


Ml.l 


11313 


N20E88 


8 


28±10 


0.45±0.29 


3.0% 


104 
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17:21 


18:08 


18:50 


5340 


Ml. 4 
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8 
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1.3% 
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M4.3 


11339 
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16 
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3.3% 
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2.1% 
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14:01 
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Ml. 9 


11389 


S25E70 


11 


33± 4 


0.43±0.05 


3.0% 


123 


2011-12-29 


21:43 


21:51 


21:59 


960 


M2.0 


11389 


S25E67 


14 


59± 


0.98±0.14 


1.6% 


124 


2011-12-30 


03:03 


03:09 


03:13 


600 


Ml. 2 


11389 


S25E67 


9 


27± 2 


0.30±0.13 


4.0% 


125 


2011-12-31 


13:09 


13:15 


13:19 


600 


M2.4 


11389 


S25E46 


13 


38± 6 


0.25±0.06 


3.6% 


126 


2011-12-31 


16:16 


16:26 


16:34 


1080 


Ml. 5 


11389 


S25E42 


15 


51± 3 


0.58zb0.08 


2.9% 


127 


2012-01-14 


13:14 


13:18 


13:20 


360 


M1.4 


11401 


N14E88 


8 


21± 5 


0.10±0.08 


3.0% 


128 


2012-01-17 


04:41 


04:53 


05:07 


1560 


Ml.O 


11401 


N18E53 


13 


38ib 2 


0.52ib0.02 


1.5% 


129 


2012-01-18 


19:04 


19:12 


19:27 


1380 


Ml. 7 


11401 


N17E32 


16 


59±11 


0.17±0.04 


2.5% 


130 


2012-01-19 


13:44 


16:05 


17:50 


14760 


M3.2 


11402 


N32E27 


38 


76±11 


0.45±0.10 


1.9% 


131 


2012-01-23 


03:38 


03:59 


04:34 


3360 


M8.7 


11402 


N33W21 


37 


111± 9 


0.98±0.09 


1.2% 


132 


2012-01-27 


17:37 


18:37 


18:56 


4740 


XI. 7 


11402 


N33W85 


40 


114± 8 


0.95±0.07 


1.7% 


133 


2012-02-06 


19:31 


20:00 


20:17 


2760 


Ml.O 


11410 


N19W62 


28 


89db 5 


0.54±0.14 


1.6% 


134 


2012-03-02 


17:29 


17:46 


18:07 


2280 


M3.3 


11429 


N19W62 


10 


36± 1 


0.75±0.03 


1.6% 


135 


2012-03-04 


10:29 


10:52 


12:16 


6420 


M2.0 


11429 


N16E65 


23 


48± 2 


0.48±0.03 


1.7% 


136 


2012-03-05 


02:30 


04:09 


04:43 


7980 


XI. 1 


11429 


N19E58 


32 


74± 5 


0.44±0.03 


1.7% 


137 


2012-03-05 


19:10 


19:16 


19:21 


660 


M2.1 


11429 


N16E45 


10 


32± 2 


0.36±0.03 


2.3% 


138 


2012-03-05 


19:27 


19:30 


19:32 


300 


Ml. 8 


11429 


N16E45 


7 


26± 5 


0.34±0.18 


3.0% 


139 


2012-03-05 


22:26 


22:34 


22:42 


960 


M1.3 


11429 


N16E43 


9 


33± 3 


0.46±0.01 


2.3% 


140 


2012-03-06 


00:22 


00:28 


00:31 


540 


M1.3 


11429 


N16E42 


8 


27± 5 


0.38di0.20 


3.1% 


141 


2012-03-06 


01:36 


01:44 


01:50 


840 


M1.2 


11429 


N16E41 


8 


28± 3 


0.40±0.04 


2.3% 
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Nr 


Obs6rva.tioii 


Start 


Peak 


End 


Duration 


GOES 


NOAA 


Heliogr. 


Length 


Diff. 


Diff. 


Fit 




date 


time 


time 


time 


T(s) 


class 


AR 


position 


L(Mm) 


coeff. K, 


index /3 


mt 


142 


2012-03-06 


04:01 


04:05 


04:08 


420 


Ml.O 


11429 


N16E39 


9 


47± 6 


0.96dz0.09 


2.0% 


143 


2012-03-06 


07:52 


07:55 


08:00 


480 


Ml.O 


11429 


N17E40 


8 


41ib 2 


0.87di0.05 


2.5% 


144 


2012-03-06 


12:23 


12:41 


12:54 


1860 


M2.1 


11429 


N21E40 


15 


55ibl0 


1.07±0.16 


1.3% 


145 


2012-03-06 


21:04 


21:11 


21:14 


600 


Ml. 3 


11429 


N16E30 


10 


38± 9 


0.59±0.11 


2.0% 


146 


2012-03-06 


22:49 


22:53 


23:11 


1320 


Ml.O 


11429 


N19E32 


11 


35± 8 


0.62±0.07 


1.8% 


147 


2012-03-07 


00:02 


00:24 


00:40 


2280 


X5.4 


11429 


N18E31 


37 


100±12 


0.60di0.01 


1.9% 


148 


2012-03-07 


01:05 


01:14 


01:23 


1080 


X1.3 


11430 


N15E26 


35 


105dz 6 


0.63dz0.03 


0.8% 


149 


2012-03-09 


03:22 


03:53 


04:18 


3360 


M6.3 


11429 


N15W03 


33 


82=b 5 


0.50ib0.01 


3.6% 


150 


2012-03-10 


17:15 


17:44 


18:30 


4500 


M8.4 


11429 


N17W24 


30 


77± 5 


0.40±0.01 


1.2% 


151 


2012-03-13 


17:12 


17:41 


18:25 


4380 


M7.9 


11429 


N17W66 


36 


83± 6 


0.73±0.01 


1.3% 


152 


2012-03-14 


15:08 


15:21 


15:36 


1680 


M2.8 


11432 


N13E05 


18 


63± 6 


0.55±0.03 


2.0% 


153 


2012-03-15 


07:23 


07:52 


08:08 


2700 


Ml. 8 


11432 


N14W03 


21 


66± 6 


0.61±0.03 


1.3% 


154 


2012-03-17 


20:32 


20:39 


20:42 


600 


M1.3 


11434 


S20W25 


10 


42± 4 


0.51±0.02 


1.1% 


155 


2012-03-23 


19:34 


19:40 


19:44 


600 


Ml.O 


11445 


S20W25 


7 


23± 2 


0.25±0.01 


2.5% 



Table 2. Statistical results of observables and model parameters of 155 analyzed flare events during 

2010/05/13-2011/03/31. 



Parameter 




Range 


Mean 


Median 


Powerlaw slope 


Flare duration (GOES rise time) T{s) = 




120-8460 


1030 ± 1050 


660 


2.17 


AIA 335 time interval Taia{s) = 




336-10992 


1806 ± 1579 


1296 




Length (flare radius) L{Mm) = 




5.5-51 


19 ±10 


16 


1.96 


Peak flux GOES Fgoes{^ m-2)= 


(1 


- 69) X 10-^ 


(4.1 ±7.8) X 10-^ 


7.8 X 10-^ 


1.92 


Total flux AIA 335 A Faia{DN/s)= 


(7.2- 


- 2462) X 10^ 


(4.4 ±4.7) X 10^ 


2.5 X 10^ 


1.34 


Mean velocity v{km/s) = 




3-103 


15± 12 


11 


2.72 


Maximum velocity Vmax{km/s) = 




8-550 


80 ±85 


57 


1.83 


Diffusion coefficient ^{km s~^/^)= 




19-139 


56 ±24 


52 




Diffusion index /3 = 




0.04-1.35 


0.53 ±0.27 


0.49 




Goodness-of-fit qfu = Ar/rmax{%)= 




0.8-4.2 


2.2 ±0.7 


2.1 
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Fig. 1. — Comparison of spatio-temporal evolution models: Logistic growth with parameters tL = 1.0, Too = 
1.0, TG = 0.1, sub-diffusion (/3 = 1/2), classical diffusion (/3 = 1), Levy flights or super-diffusion (/3 = 3/2), 
and linear expansion (r oc t). All three curves intersect at t = ti, and have the same speed v = {dr/dt) at 
the intersection point at time t = tL- 
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Fig. 2.— The flare event #28 was observed on 2011 Mar 7, 19:43-20:58 UT, with AIA/SDO 335 A, with 
GOES time profiles (top panel), the total flux (second panel), the spatio-temporal evolution of the radius 
r{t) = y^A{t)/7T of the time-integrated flare area A{t) for two thresholds, Fthresh = 100,200 DN/s (third 
panel; histogrammed), fitted with the anomalous diffusion model (third panel; solid curve), and 5 snapshots 
of the baseline-subtracted flux (fourth row) and highpass-filtered flux (bottom row) , with the threshold flux 
Fthresh = 100 DN/s showu as contour. 
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Fig. 3. — A selection of analyzed flare events, showing the GOES 1-8 A time proflle (dotted curve), the AIA 
335 A total flux (dashed curve), the measured radius of the flare area for two flux thresholds (histograms), 
the best-fit diffusion model (solid curve), and prefiare-subtracted AIA 335 image at the end of the flare, with 
the contours of the flux threshold of Fthresh = 100 DN/s indicated (right panels). 
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Fig. 4. — Distributions of various observed parameters (L, T, Fg, v, Vm, Fa) in log-log format with powerlaw 
fit (left and middle column), and best-fit model parameters {k, (3, Ar/rm) with linear histogram and mean 
and standard deviation indicated (right column). 
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Fig. 5. — Correlations between observed and best-fit parameters. The linear regression fits are shown for 
y{x) and x{y) (thin lines), and for the arithmetic mean (thick lines). 



